Effect of nearest- and next-nearest neighbor interactions on the spin-wave velocity 
of one-dimensional quarter-filled spin-density-wave conductors 
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We study spin fluctuations in quarter-filled one-dimensional spin-density-wave systems in pres- 
ence of short-range Coulomb interactions. By applying a path integral method, the spin-wave 
velocity is calculated as a function of on-site (U), nearest (V) and next-nearest (V2) neighbor-site 
interactions. With increasing V or V2, the pure spin-density- wave state evolves into a state with 
coexisting spin- and charge-density waves. The spin-wave velocity is reduced when several density 
waves coexist in the ground state, and may even vanish at large V . The effect of dimerization along 
the chain is also considered. 
PACS Numbers: 72.15.Nj, 75.30.Fv 



I. INTRODUCTION 

Organic conductors of the tetramethyltetraselena- 
fulvalene (TMTSF) and tetramethyltetrathiafulvalenc 
(TMTTF) salts family often exhibit density-wave (DW) 
instability at low temperatureJina Recent experiments 
have shown that a 2k F spin-density wave (SDW) may 
coexist with a Ak F and/or a 2k F charge-density wave 
(CDW)rrtl (The quantity k F denotes the one-dimensional 
Fermi wave vector and 2k F is the nesting wave vector for 
the SDW.) Furthermore, these CDW's seem to be of pure 
electronic origin without any (significant) contribution 
from the lattice. 

This unusual ground-state can be understood on the 
basis of a mean-field theory for a quarter-filled one- 
dimensional system in the presence of several kinds of 
Coulomb interaction. Within an extended Hubbard 
model with on-site (U) and nearest-neighbor (V) inter- 
actions, it has been shown that a 4k F CDW may cct 
exist with the 2k F SDW when V is strong enough.tl 
When the next-nearest-neighbor interaction (V2) is also 
taken intcuaccount, three different ground states can be 
stabilizedfTO (i) a pure 2k F SDW at small V and V 2 , 
(ii) coexisting 2fc F SDW and 4k F CDW at large V, (iii) 
coexisting 2k F SDW, 2k F CDW and 4fc F SDW at large 
V2. Although the SDW instability is driven by the on- 
site repulsive interaction U, the nearest and next-nearest 
neighbor interactions play a crucial role for the appear- 
ance Of CDW'S. r-| |— - 

Following the standard analysis ,E!nLZl fluctuations 
around the mean-field ground-state have been studied. 
For a quarter-filled system, commensurability effects with 
the underlying crystal lattice pin the DW's and produce 
a gap in the sliding modes.E3 Surprisingly, this gap van- 
ishes at the boundary between the pure 2k F SDW and 



the coexisting 2k F SDW and 4fc F CDW.E1 The spin- 
wave modes have been, studied only within the Hubbard 
model (V = V% = 0)£3 The spin- wave velocity decreases 
monotonically with increasing U, in qualitative agree- 
ment with the. exact solution of the one-dimensional Hub- 
bard modelEil 

In this paper, we study the spin- wave modes in pres- 
ence of the nearest and next-nearest neighbor interac- 
tions (V, V% ^ 0) . We consider a one-dimensional system, 
assuming that long-range order is stabilized by (weak) in- 
terchain coupling. Gar ^analysis is based on a functional 
integral formulationESca which allows a simple treat- 
ment of the spin- wave modes even in the presence of these 
interactions. The electron-electron interaction is treated 
within (Hartree-Fock) mean-field theory, while the SU(2) 
spin rotation symmetry is maintained by introducing a 
fluctuating spin-quantization axis in the functional in- 
tegral. Transverse spin-wave modes then correspond 
to fluctuations of the spin-quantization axis around its 
mean-field value. 

In Sees. II and III, we extend the derivation of Ref. |25| 
from the incommensurate, to the commensurate case. 
Previous mean-field resultsa are recovered within a saddle 
point approximation. Then we derive the effective action 
of the spin-wave modes and obtain the spin-wave veloc- 
ity. In Sec. IV, the spin-wave velocity is calculated as a 
function of V, V% and the dimerization along the chain. 
Section V is devoted to discussion. 



II. PATH INTEGRAL FORMULATION 



We consider a one-dimensional electron system at 
quarter-filling with dimerization along the chain. Within 
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the extended Hubbard model, the Hamiltonian is given 
by 



H = Hq + Hj , 



(2.1) 
(2.2) 



<7,n,n' 



Hi = U ^n nT n„; + V ^(V^nXV'j.+iV'n+i) 

n n 

n+2j 



C7 



(2.3) 



where ip n = (ipnt, ^nj.)*, ^ncr = V^Vw, and V^L is tnc 
creation operator of an electron with spin cr(=t, |) at the 
lattice site n. The transfer integral in the kinetic term 
Hq is defined by 

( t - for n' = n + l , 

W = < t + (-l) n t d for n' = n - 1 , (2.4) 
I otherwise , 

where a finite is due to the dimerization. The inter- 
action Hamiltonian is expressed in terms of the Hubbard 
interaction U and the density-density interaction V nn i 
defined by 



U/4 for n' = n , 

_ J V/2 for ri = n ± 1 , 

Va/2 for n' = n±2 , 

otherwise , 



(2.5) 



where V (V2) is the coupling constant for nearest (next- 
nearest) neighbor-site interaction (U, V, V2 > 0). 

In order to derive the effective action for the spin- wave 
modes, we write the partition function Z as a path inte- 
gral: 



Z = / V^Vi\) e 



5 



r/r 



(2.6) 



(2.7) 



where the action S is a function of the Grassmann vari- 
able tjj. t is a Matsubara time- varying between and 
1/T. Following Refs. ^3] and |25|, we now introduce the 
new field <fi defined by 



Ipn = Rn4>n 

R n a z Kl = a- ■ n„ 



(2.8) 



where R n is an SU(2)/U(1) unitary matrix and n„ is 
a unit vector which gives the direction of the spin- 
quantization axis at site n and time r for the field </>. 
Substituting Eq. ( |2.8D into Eq. (2.7), the action is rewrit- 
ten as S = So + Si, where 



S 



Sr 



j dr^^Udr-^ + RldrRn 

- i'nRntnn' Rn' 4>; 
n,n' 

j dT ) ~J zZ An + Z-J P^Vnn'Pc 
\ n n,n' 



(2.9) 
(2.10) 



Pen = ^n^n and p sn — (fyo z (f> n are the charge- and spin- 
density operators. The quantities cr x , a y and er 2 are Pauli 
matrices. Note that Si is invariant under the transfor- 
mation if) — > <j), since the interaction is invariant with 
respect to spin rotations. It is convenient to rewrite the 
action as 



S = 



dT I zZ ^n(^T - M - Aon)4>n 
I n 

^2 0nW ex P _i j 



dlA x i 4> n 



U 



TzZ^n + Yl P™Vn 



(2.11) 



where the SU(2) gauge fields Ao and A x arc defined by 

A Qn = -R n d T R n , (2.12a) 
R n R n+s , (S = ±1) . (2.12b) 



pn+8 

expl -1 j dlA x i 



J 



The lattice spacing is taken— .as unity. Using the 
Stratonovich-Hubbard identity^, the interaction part of 
the action is rewritten as (note that U, V, V2 > 0) 



expl ~22 drpcnVnn'pcn' = / £>A C exp ~^Z A cn V n ^, A cn , + 2iJ^ / drA 

\ n,n' J \ n ' n> n 

ex p(jE / dT p 2 ^ = J pa * cx p(-tjE / drA2 sn+ }Z J drA 



cnHcn I i 



sn f J sn I i 



(2.13) 
(2.14) 
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where A c n an d A srl are (real) auxiliary fields. By us- 
ing Eqs. ( 2.13 ) and ( 2.14 ), the final form of the partition A c „, we obtain the self-consistent mean- field equations 
function is given by 



Z = J VA C VA S J Vn J V$V<t> e - {So+Sl) 
S = J dr < <f>i,(d T - (Jt- -4on)0n 



(2.15) 



(Psn )mF ' 



(3.1) 
(3.2) 



E <Pitnn> exp 
" 1 



dlA x l 4> n > 



(2.16) 



where V - ! = FT 1 . 



jj ^sn ^snPsn 2z A cn p c 



(2.17) 



By minimizing the free energy with respect to A sra and 
i se 

_ U 
sn ~ 2 

The average ( )mf is to be calculated with the mean-field 
action 

n n,n' 

+ J dr |^ ^ (9 T - /i - 2iA cn - A sn cr z ) <^„ 



E 0nW<£n' 



At quarter-filling, the mean-fields (p. 



(3.3) 



III. EFFECTIVE ACTION FOR THE SPIN- WAVE 
MODE 

In this section, we derive the action corresponding to 
the spin-wave modes at quarter-filling. First, we repro- 
duce the mean-field result of Ref. || within a saddle-point 
approximation. Then we consider transverse spin fluctu- 
ations arising from the dynamics of the spin-quantization 
axis. 



Sn /MF & n d 

(Pc?t.)mf are periodic with a periodicity of four lattice 
spacings. They can be written as 



(A 



sn/MF 



E Sv 



m— 
3 



(Pen) MF — 

m=0 



where Qo = 2kp = it/2. Since (p c 



imQon 



o irnQ n 



(3.4) 
(3.5) 



and (p s 



A. Mean-field solution 



'cii/mf auu \rsn/MF 

real quantities, one finds D = Dq,Di = D\,D 2 = D 2 
and S = S$,Si = SS,5 2 = S|. In Eqs. Q and 
d3.5[ ), S'o = due to the absence of ferromagnetism an d 
The standard mean-field solution is recovered from a Dp = 1/2 for a quarter-filled band. From Eqs. (3.1)- 



saddle-point approximation with n = z at each lattice 
site. One then has R n — 1 and A$ — A x = 0. 



5h, the final form of the mean-field action is obtained 



J 



Smf = {3N 



~ fd^il 2 - \Si\ 2 ) S\) V(\ Dl) V 2 {\ 2\D^ + Dl) 



/^{E4 ( d r 



/i + — + V + V 2 - 2t cos k 



E 4 [j(Di - Si*,) - 2V 2 Di) 4> h - Qo + c, 
E 4 G°2 - S 2 *z) - 2VD 2 + 2V 2 D 2 - 2ii d sin fc) <fe_ 2Qo | 



(3.6) 



where <^>ft = (l/VN)^2 n e~ lkn <t>n an d AT is the num- 
ber of lattice sites. The action ( |3.6| ) agrees with the 
mean-field Hamiltanian obtained previously by the con- 
ventional methodic 
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B. Fluctuations 



In the long- wavelength limit, collective modes can 
be separated into sliding (charge) modes and spin-wave 
modes. In this paper, we consider only transverse (acous- 
tic) spin- wave modes (i.e. magnons). These modes show 
up in the fluctuations of the unit vector field n. They 
do not couple to charge modes and gapped spin-wave 
modes. We shall make the following two approximations: 
(i) We neglect the coupling to long-wavelength spin fluc- 
tuations [A s (q) with \q\ <C Q ]. In the Hubbard model 
{V = V2 = 0), this coupling is known to renormalize the 
spin- wave velocity Jay the factor [1 - UN(0)} 1/2 in the 
weak-coupling limito [N(0) is the density of states at the 
Fermi level] ; (ii) We also neglect any possible coupling to 
spin fluctuations, at wave- vector 2Q + q [A S (2Q + q) 
with \q\ < Qo}M 

When two SDWs coexist in the ground-state, our for- 
malism can only yield the "in-phase" modes where the 
two spin-density waves oscillate in phase. It misses the 
modes where the oscillations are out-of-phase.Ej These 
modes are gapped and do not couple to the "in-phase" 
modes considered in this paper. 

Before proceeding with the spin-wave mode analysis, 
let us discuss the limit of validity of our approach. The 
spin-wave modes will be obtained by expanding about 
the (Hartree-Fock) mean-field state. Such an approach 
should hold (at least qualitatively) as long as the inter- 
action is smaller than the bandwidth, i.e. U, V, V2 < 
At. Nevertheless, it does not necessary break down in 
the strong-coupling limit. In the context of the two- 
dimensional Hubbard model, Schrieffer et al. have shown 
that an RPA analysis of the fluctuations about the mean- 
field state in the limit U 3> t agrees with the conclusions 
obtained from thp-Heiscnberg model with exchange con- 
stant J = u 2 /uH 



J 



Another limitation of our approach comes from the 
analysis of the fluctuations of the unit vector n. As 
will become clear below, the main assumption is that 
n is a slowly varying field, thus allowing a gradient ex- 
pansion. Whereas this assumption is perfectly valid in 
the weak-coupling limit (U,V,V2 ^ At), it breaks down 
in the strong-coupling limit. In the latter, one should 
write n„ = n^ low + cos(n7r/2)L n where n^ low is a slowly 
varying field and L„ a small perpendiciilap-|Component 
(L„ • < low = and |L„| < \nf° w \ ~ 1)BE3 The effec- 
tive action of the spin- wave modes, 5 , e ff[n slow ], is then 
obtained by integrating out both the fermions and the 
(small) transverse component L„. For V = Vi = 0, 
this allows to interpolate smoothly between the weak- 
coupling regime and the strong-couplingr*egime which is 
well described by the Heisenberg modelEa 

Long-wavelength transverse spin fluctuations corre- 
spond to fluctu ations of the SU(2) gauge fields Aq and 
A x [Eqs. fl2.12a| ) and ( [2.12b| )l which are rewritten as 



Aon — ^2 ^0n°V ) 
A X n = ^2 ■ 



(3.7) 
(3.8) 



From eqs. (2.11), (p\q), (3.7) and (3.8), we write the ac 
tion of the spin degrees of freedom as 

S = Smf - E J dT(t>nA 0n (j) n 



t nn , exp —i I dlA x 



(3.9) 



To order 0(A 2 ) we obtain 



S = Smf - E J dr^Aon^n - E J ^nn'^n ^ ^ n ~ ^ xn + ^xn') ~ 2^" 1 ) 



Si 



MF 



EE E 

n fj,—0,x i/—x,y,z 



dTj»A»+s: 



dia 
x 1 



(3.10) 



where j£ n , j% n and a are given by 



JOr, 



5=±l 



5 " ia = 5EX>«»' / dT <tiwt>n>Ai n Ai n 



(3.11) 
(3.12) 

(3.13) 



The second term of Eq. (3.10) denotes the coupling of 
the gauge field A" to the spin current (j% n ) and spin 



density (jftj. The last toon of Eq. ( pU0| ), Sf a , is the 
diamagnetic contributions 

The effective action of the gauge field is obtained by 
integrating out th e ferm ions in the part ition function. By 
substituting Eq. fl3.10| ) into Eq. (§T5|), one obtains the 
effective action up to 0(A 2 ) as 

S eS [A»] = (S^) MF - E / dT(j; n ) MF A» n 
"^E E Jdrdr'A^r) 

n n 1 11 11 1 u 1/' 



71,71 {£,{£' ,V.V 



A 



xn.^ ; (n,T,n',T')^(r') , (3.14) 



where 



J u, J u'(=0,x) 



(3.17) 



n 3 ^(n,r,n / y) = (^(r)^ n/ (0) , (3.15) 
(5f) MF = J E *»«' / dr(4^> MF K„) 2 • (3-16) 



The quantity II .„ .„/ denotes the current-current corre- 

lation function in the mean-field state. We, note that 
(j^ n ) MF = in the long- wavelength limited and that 



(K) 



MF 



(3.18) 



MF 



where (K) MF is the mean value of the kinetic energy 
per site in the mean-field state, q = (q, iCl) and f2 is 
a bosonic Matsubara frequency. The quantity (q) 

is the Fourier t ransf orm of Eq. ( 3.1 5| ) with respect to n 
and r. In Eq. (3.17), it can be evaluated at q = since 



A v is of the order O(V). To order 0(V 2 )), we obtain A £ « °( V )' Note that = 11 



n ./..'./., 







->cff 



|£{<*>mf E l^(9)l 8 

5 I v=x,y,z 



and IT-v ,v (a) _ „ = 0. Taking the continuum limit 
n — > £ (with £ a real continuous variable) and writing 
A^ n — Ap(£, r), the effective action ( ft.!7| ) is rewritten as 



<5 C ff 



-5E E {wmfK®i 2 +Ei^®I 2ii «} 

{ i/=x,y u=x,y 

- \ E + 'I,, ) + n«l^(9)l 2 } , 



(3.19) 



IljHjy . Here 



where Ui*^ = Tlju^fq — 0) and IL^ 
we note the identities (K) MF + = and n,-*^ = 

whic h ca n be deduced from the gauge invariance of 
Eq. (|]^) (Appendix A). We have verified numerically 
the validity of these identities. Finally, noting thatEJ 



■=x,y 

^ 2 (£,r) = -i(d r n) 2 , 



(3.20) 
(3.21) 



v—x^y 



we obtain the following final - 
action of the spin-wave modes 



ptntession for the effective 
ESELl (Appendix A): 



S cS = \ J d£dr { X (d T n) 2 + P (d^) 2 } , (3.22) 

where x ancl P are the uniform transverse spin suscepti- 
bility and the spin stiffness, respectively: 



_ 1 q q \MF 



1 . 



q =a ^i'oi'a > 



{v = x,y) 



P = -4((^>MF + n «) » (v = x,y) . 
From Eq. ( |3.22| ) we deduce the spin- wave velocity 



(3.23) 
(3.24) 



(3.25) 



In the incommensurate case, —i-0 in the weak cou- 

pling limit so that p = — {K) MF /4£3 As shown in the 
next section, gives rise to a contribution of the 

same order as (K) MF in the quarter-filled case when the 
on-site i ntera ction U is of the order of the bandwidth. 
In Eqs. d3~23| ) and ( gjgg ), j£ and K can be expressed as 
[<pk = (<t>kh falY, v = x,y\ 



jo (?= 0) = -^E^!^ 



ix(? = °) = ^E (2* sin. fc^ov^k 

- 2ii d COS fc ^0- v ^fc+2Qa J 

^ = ^E(- 2<cosfc ^ 



2^ d sinfc0^ fc+ 2Q t 



(3.26) 



(3.27) 



(3.28) 
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IV. SPIN- WAVE VELOCITY 



In this section, we evaluate the spin-wave velocity at 
zero temperature (T = 0). We take t = 1 and calculate 
the velocity normalized to its value at V = V2 = and 
i d = 0. 

The phase diagram of the present model as a function 
of V and V2 is shown in Fig. 1 for U = 4 and id = (solid 
curve). El For small V and V2, there is a pure 2/cf SDW 
state (region I). A large V induces a phase with both a 
2fcp SDW and 4fcp CDW (region II), while in the presence 
of a large V2 there is coexistence between a 2fcp SDW, 
a 2fc F CDW and a 4fc F SDW (region III). The dashed 
curve denotes the boundary at which a first order tran- 
sition occurs between II and III. The dash-dotted curve 
shows the phase diagram for t& = 0.1. The sliding modes 
are gapped in all three regions. However, the charge fluc- 
tuations become gapless at the transition between I an d 
II. We discuss below the spin-wave velocity [Eq. (||5|)] 
as a function of V and V2 for both td = and id 7^ 0. 



A. U dependence (V = V2 = and t d — 0) 



The spin stiffness p and the susceptibility x are shown 
in Fig. 2(a) as a function of U for V = V2 = and 
id = 0. Both p and x are almost constant for small U 
and decrease monotonically for large U. The inset shows 
the corresponding U -de pend ence for (K} MF and Hj*jm 



0.4 



which determine p [Eq. (3.24)]. A behavior similar to 
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FIG. 1. Phase diagram in the V — V2 plane for U = 4 
and td = (Ref. |^). The three different regions corre- 
spond to: (I) pure 2k F SDW, (II) coexisting 2k F SDW 
and 4fc F CDW, (III) coexisting 2fc P SDW, 2fc P CDW 
and 4fcp SDW. The dash-dotted curves denote the corre- 
sponding boundaries for id = 0.1. 
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U 

FIG. 2. (a) U dependence of p and x f° r V = F2 = 
and td = 0. The inset shows the U dependence of (K) MF 
and Iljgji . (b) U dependence of v (solid curve) and 
v[l — 2\U] 1 ^ 2 (dashed curve). The open circles denote 
the exact result for the one-dimensional Hubbard model 
(Ref. §). 

the incommensurate case is seen for U < 2: Hp.,* is van- 
ishingly small, and x, p and (K) M -p are almost constant 
with respect to U. The limiting values for small U are 
given by x = l/(2%/2vr) ~ 0.113, (K) MF = -2^/^ ~ 
—0.90, and v = y/2. The variation of these quantities 
for U > 2 comes from the effect of commensurability at 
quarter-filling. 

In Fig. 2(b) (sol id curve), we show the spin-wave ve- 
locity v [Eq. ( 3.25| )], which is almost independent of U 
although slightly suppressed at large U. Here we note 
that we have neglected the coupling to long-wavelength 
spin fluctuations. In the Hubbard model (V = V2 = 0), 



G 



the spin- wave velocity v = (p/x) 1 ^ 2 becomes (p/x) 1 / 2 [l— 
2xU] x ' 2 when this coupling is taken into account within 
the RPA.EZl One obtains 1 - 2 X U = 1 - UN[0) in the 
weak-coupling limit where N(0) = l/^/2n at quarter- 
filling. In Fig. 2(b), we show v and v[l — ZxU] 1 / 2 (dashed 
curve). The open circles denote the exact result for the 
one-dimensional Hubbard modelEd For [/ < 2, the RPA 
result turns out to be a good approximation, while the 
difference becomes noticeable at larger U. Nevertheless 
we use v = (p/x) 1 ^ 2 as a nrs t step to examine the spin- 
wave velocity as the function of V and V 2 , The present 
calculation is performed by choosing [7 = 4, which leads 
ton - 1.29 for V = V 2 = and t d = 0. 



B. V dependence (V 2 = and td = 0) 



Now we consider the V dependence of the spin-wave 
velocity for V 2 — 0, t d = and U/t = 4. Contrary 
to the weak-coupling limit which can, be studied analyt- 
ically as the incommensurate case£3 this intermediate 
coupling regime requires numerical calculation. Figure 3 
shows the ^-dependence of v, x an d P ( au quantities are 
normalized to their value at V = V2 = and t<j = 0). 
The arrow indicates the critical value V c = 0.34 separat- 
ing regions I (Si 7^ 0) and II {S\ 1 D 2 7^ 0). In region 
II (V > V c ), both p and \ decrease for decreasing V. 
The stronger decrease of p results in a decrease of the 
spin-wave velocity. For large V, both the spin stiffness 
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V 

FIG. 3. V dependence of the spin- wave velocity ve- 
locity v, susceptibility \ and spin stiffness p for (7 = 4, 
t d = and V 2 = 0. v°(= 1.286), x°{= 0.103) and 
p°(= 0.171) are the values for V = V 2 = 0. There is a 
small jump at the critical value V c = 0.34, which is shown 
by the arrow. The inset shows the F-dependence of the 
order parameters Si (2fc F SDW) and D 2 (4fc P CDW). 



T 1 1 J 1 1 1 J 1 1 1 1 1 1 r 




FIG. 4. V 2 dependence of v, x an d p for U = 4, 
td = and V — 0. There is a cusp at the critical value 
V2c ~ 1-32 corresponding to the transition from state I to 
state III. The inset shows the ^-dependence of the order 
parameters Si, £>i (2fc F CDW) and S 2 (4fc F SDW). 

and the spin-wave velocity vanish. It seems that the de- 
crease of v in region II mainly comes from the reduction 
of kinetic energy due to the formation of the 4fcp-CDW. 
Note that the spin-wave velocity is discontinuous at the 
critical value V = V c . The small jump at V c originates in 
the discontinuity of Si axid D 2 (see inset of Fig. 3) which 
is found only for id = 0.EI 



C. V 2 dependence (td = 0) 

In this section, we analyze the V 2 dependence of the 
spin-wave velocity for U = 4, t d = and different values 
of V. 

Figure 4 shows v/v°, x/x° an d p/p° in the case V = 
(the inset shows Si, D\ and S 2 as a function of V 2 ). 
There is a transition between regions I and III at the 
critical value V 2c . v/v°, x/x° an d p/p° are constant for 
V 2 < V 2c , and decrease for V 2 > V 2c (note that v ac- 
tually slightly increases at large V 2 ). However, all these 
quantities remain finite in the limit of large V 2 . This is 
to be contrasted to the large-]/ limit (region II) where 
the spin- wave velocity vanishes (Fig. 3). Such a behav- 
ior can be understood as follows. For V2 — ► 00 (region 
III), the spin- and charge-density waves in the ground- 
state are of the type (T, 1,0,0) and (1,1,0,0), respectively. 
Our numerical calculation shows that this behavior al- 
ready shows up for V 2 /t ~ 4. In this limit (V 2 /t > 4), 
the one-dimensional chain divides into independent two- 
site clusters. For this problem, one can find the exact 
expression of the spin- wave velocity (Appendix B): 
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(inset) exhibit a similar behavior (Fig. 5(a)). They are 
constant in region I, and increase (decrease) in II (III) 
when V2 increases. Figure 5(b) shows the spin-wave ve- 
locity v/v° which turns out to be mainly determined by 
p/p°. Except in region I and for large values of V27 v 
varies strongly as a function of V2 ■ 

Here we comment on the fact that v remains finite 
at large V%. Within the mean-field treatment, which is 
expected to be valid for a moderate coupling between 
chains, both \ an d P remain finite at large V2. On the 
other hand, for one-dimensional systems it is known from 
bosonization that \ vanishes at large V2 due to the for- 
mation of a spin gapEJ Thus, we expect our mean-field 
analysis in region III of Fig. 1 to break down when the 
interchain coupling becomes sufficiently small. 



0.5 



1 1 1 1 1 1 1 1 1 
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/ 2 /v 








/ V=3 




- s 


(b) ; 
. 1 . . . 



2 

v 2 



FIG. 5. (a) V 2 dependence of x for V=0 (Fig. 4), 1, 
2 and 3. There is a jump due to the first order transi- 
tion between II and III indicated by the dashed curve 
in Fig. 1. The inset shows the ^-dependence of p. 
(b) Fa-dependence of v for V=Q (Fig. 4), 1, 2 and 3. 
For V=l, there is a small jump at Fj=0.66, which corre- 
sponds to the transition between II and I. This disconti- 
nuity can also be seen in Fig. 3. 



1/2 



/v° = (t- t d )/v° 



(4.1) 



For U = 4, v° = 1.286, so that v/v° = 0.777. For 
Va/i = 4 and 13, the numerical calculation gives v/v° — 
0.763 and 0.776, respectively, in excellent agreement with 
the analytical result of the two-site problem. 

Now we consider the V"2-dependence of x/x°> p/p° an d 
v/v° for V = 0, 1, 2 and 3. For V = 1, there is first a 
transition from region II to region I, and then a transition 
from I to III. For V = 2 or 3, there is a single transition 
occuring between II and III. The ratio x/x° and p/p° 
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FIG. 6. (a) V dependence of v for U = 4, V 2 = 
and td=0, 0.1, 0.3 and 0.5. (b) V2 dependence of v for 
U = 4, V = and i d =0, 0.1, 0.3 and 0.5. 
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D. Effect of dimerization 

Finally, we consider the effect of dimerization on 
the spin-wave velocity v. Figure 6(a) shows the In- 
dependence for U = 4, V2 = and t^ = (solid curve), 
0.1 (dotted curve), 0.3 (dashed curve) and 0.5 (dash- 
dotted curve). The effect of dimerization is large in re- 
gion I, but rather small in region II. A finite increases 
the band gap. This induces a suppression of Hj*j* and 
p, and leads ultimately to a reduction of the spin-wave 
velocity. We note that the reduction of Si and D 2 in re- 
gion II by the dimerization has little effect on v, since the 
dependence of Si and D2 on dimerization is very small 
for V < A. 

Figure 6(b) shows the ^-dependence of v for U = 4, 
V = and id — (solid curve), 0.1 (dotted curve), 0.3 
(dashed curve) and 0.5 (dash-dotted curve). The effect 
of dimerization is noticeable in both regions I an d III . 
The limiting behavior for large V 2 is given by Eq. (4.1). 
In that limit, the SDW exists for U > 2{t — t<j) and the 
spin- wave velocity v depends only on t and t^. 



V. CONCLUSION 

In conclusion, the nearest and next-nearest neighbor 
interactions strongly affects the spin-wave velocity in the 
intermediate coupling regime U ~ At. Our main results 
are as follows, (i) In the pure SDW state (region I), the 
spin-wave velocity v is independent of the nearest (V) 
and next-nearest (V2) interaction (Fig. 3). (ii) For coex- 
isting 2kp SDW and Akp CDW (region II), v decreases 
(increases) as a function of V (V2) [Figs. 3 and 5(b)]. It 
is slightly discontinuous at the transition between I and 
II and vanishes (as well as the spin stiffness) at large V 
(Fig. 3). (iii) For coexisting 2k F SDW, 2fc F CDW and 
4fcp SDW (region III), v is suppressed by V^- It tends 
to a finite value at large V2 [Figs. 4 and 5(b)]. (iv) The 
dimerization decreases the spin-wave velocity [Figs. 6(a) 
and (b)]. 



J 



As discussed in Sec. III.B, our approach is limited to 
the weak to intermediate coupling regime and should hold 
when U, V, V 2 ^ At. In the half-filled Hubbard model, a 
strong coupling is known to reduce the spin- wave velocity 
from v = 0(t) tov= O(J) (with J = At 2 /U < t). We 
also expect a decrease of the spin-wave velocity in the 
more general case we have studied when U, V, V 2 become 
larger than At. Therefore, our main conclusion (a reduc- 
tion of the spin- wave velocity by the interactions V, V2) is 
likely to be strengthened by strong coupling effects. The 
Stoner factor (l—2xU) 1 ^ 2 , which arises from the coupling 
to long-wavelength spin fluctuations, was not considered 
in our analysis. It leads to a decrease of v when the on- 
site interaction U increases. Whether the Stoner factor 
depends on the interactions V and V2 remains an open 
question. 

In the cCkHipounds that have been studied 
experimentallyJJ'Q the Bechgaard salts (TMTSF) 2 PF 6 
and (TMTSF) 2 AsF 6 , and the Fabre salt (TMTTF) 2 Br, 
the electron-electron interaction is expected to be in the 
intermediate coupling regime (U ~ At). Furthermore, es=. 
timates by MilatJ and quantum-chemistry calculationsEII 
have revealed the finite-range part of the Coulomb poten- 
tial, the first-neighbor interaction V being equal or even 
larger than U/2. We therefore think that our conclu- 
sions are relevant to the Bechgaard-Fabre salts studied 
111 Refs. I and |. 



ACKNOWLEDGMENTS 



We (Y.T. and Y.S.) thank H. Sakanaka for useful dis- 
cussion. This work was financially supported by Uni- 
versite Paris-Sud, France and a Grant-in- Aid for Scien- 
tific Research from the Ministry of Education, Science, 
Sports and Culture (Grant No. 09640429), Japan. It was 
also supported by Core Research for Evolutional Science 
and Technology (CREST), Japan Science and Technol- 
ogy Corporation (JST). 



APPENDIX A: DERIVATION OF EQ. (3.22) 

We rewrite Eq. ( 3.19[ ) as 

S cS = \J d^dr { X (d T nf +p(^n) 2 } - \ ^ { ((K) MF + \A%{q)\* + Il jSjS \A^q)\ 2 } 

q 

where p and x are given by 

P= -\(( k )mf + u j i ) = -\(( k )mf + ) » 

and 



(Al) 



(A2) 



X — (S V S V ) MF — -^N-jgjv , 



v = x,y 



(A3) 
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To show that the second term of Eq. (Al) vanishes, we use the invariance of the action under the gauge transformation 
^-u(C) T ) — * ^u(£i r) + i<9 M A(£, t) ([i = £ or t)Ej. This transformation corresponds to a rotation of n MF = z around 
the z axis and does not change the state of the system. The invariance of the action in this gauge transformation 
implies 



JE{(Wmf+IW 



~q 2 x \A(q)\ 2 -iA*(q)q x A(-q) 



n 



±n 2 \A(q)\ 2 + iA*(q)nM-<i) 



= o 



(A4) 



Since Eq. (A4) should be valid for an arbitrary function 
A, we deduce 



(K) 



(A5) 
(A6) 



which lead to the vanishing of the second line of Eq. (Al). 

Equations (A5) and (A£) can also be obtained from 
the U(l) electromagnetic field gau ge in variance. Noting 
that Iljz jz = Hjo jo , Eqs. ( A5) and ( A6) can be rewritten 
as 



( k )mf + n jSiS = ' 
. 



-'O-'O 



(A7) 
(A8) 



We recognize here the components of the polarization 
tensor for the usual U(l ) electromagnetic gauge field. 
Equations (A7) and (A8) follow from (electromagnetic) 



gauge invariance 



APPENDIX B: LIMITING CASE OF LARGE V 2 



When V2 — ► oo, the mean- field solution in region III 
corresponds to that of a half-filled two-site system given 
by 



H = -(t- i d ) E ( C l C ^ + H - c - 

+ U (ni^n ll +n 2 ^n2i) + V^ni a n 2 a' , (Bl) 



where n 1<y (n 2<7 ) = C\ a C la [C\ a C 2a ) and C\ a (C\ a ) dc- 
note the creation operators of an electron at site 1 (2) 
with spin a. The mean-field equations are given by 



(B2) 
(B3) 
(B4) 



where the average ( ) is performed with the mean-field 
Hamiltonian 



(J (7 

E(c i t CT C lCT )sgn(a) = A , 
WcLC 2CT )sgn(a) = -A 



%f = -(t - t d ) (C'lC 2 a + h.c^j 

+ E 



^ + V~sgn(a)^A)ciC 1<7 



| +V + sgn(a)jA)ctC 2o 



(B5) 



From Eqs. (p32|), (B2), (B4) and (B5), the self-consistency 
equation for A is expressed as 



1 



U/2 



y/[{U/2)A]* + {t-t d y 



(B6) 



where a = U/2 + V at half-filling. The solution of 
Eq. ( |B6| ) is obtained as 



±\ 1 



2{t-t d ) 
U 



(B7) 



for U/(t — td) > 2. By using Eq. (B7), we compute the 
uniform transverse spin susceptibility {\') and the spin 
stiffness (//): 



1 



• E 

n,n' — 1,2 

2U 



4 

2{t-t d ) 
U 3 



Oo»j o >0) L=o 



'K') + n' 



"-'■""is? 



2{t-t d f 
U 3 



(B8) 



(B9) 



where the kinetic energy per site ((K'}) and the spin 
current-current correlation function (n^-x) are given by 



2{t-t d f 
U 



\ E mn)f x {n')) | in=0 



2 

8ft - t d ) 
U 3 



n,n' = l,2 
4 



(BIO) 



(Bll) 



with 
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f x (n) - 



(B12) 



^EK C ^-Ct^-.) ■ (B13) 



By noting that \ = x'/2 and p = p'/2, we obtain 
the s pin- wave velocity of the one-dimensional system 
[Eq. (pjj)] in the limit V2 —* 00 as 



1/2 



(B14) 
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